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Abstract 

The quasinormal modes (QNMs) associated with the decay of Dirac field perturbation around a 
Schwarzschild black hole is investigated by using continued fraction and Hill-determinant approaches. 
It is shown that the fundamental quasinormal frequencies become evenly spaced for large angular 
quantum number and the spacing is given by wa+i — = 0.38490 — O.OOOOOi. The angular quan- 
tum number has the surprising effect of increasing real part of the quasinormal frequencies, but it 
almost does not affect imaginary part, especially for low overtones. In addition, the quasinormal fre- 
quencies also become evenly spaced for large overtone number and the spacing for imaginary part is 
Im(tu n+ i) — Im{uj n ) ~ —i/4M which is same as that of the scalar, electromagnetic, and gravitational 
perturbations. 
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I. INTRODUCTION 



Motivated by the beliefs that one can directly identify a black hole existence by comparing 

QNMs with the gravitational waves observed in the universe, and the study of QNMs may lead 

to a deeper understanding of the thermodynamic properties of black holes in loop quantum 

gravity 111 |2j , as well as the QNMs of anti-de Sitter black holes have a direction interpretation 

flflfl 

in terms of the dual conformal field theory [3( |i| [5(, the investigation of QNMs in black hole 
spacetimes becomes appealing recent years. Since Regge and Wheeler (J first presented the 
idea of the QNMs and Chandrasekhar and Detweller 0] worked out the numerically calculation 
of the QNMs for the Schwarzschild black hole, a great deal of effort has been contributed to 
compute the QNMs of the black holes nj- 'M\ for fields of integer spin, such as the scalar, 



electromagnetic and gravitational perturbations. 

On the contrary, the study of the QNMs of the Dirac field is very limited [2j| [2i| [2j|. Cho |2J 
studied the QNMs of the Dirac field of the Schwarzschild black hole with the help of the third- 
order WKB method, and Zhidenko \2m extended the study to the case of the Schwarzschild-de 

n 

Sitter black hole by using the sixth-order WKB method. We .27] investigated the Dirac QNMs 
of the Reissner-Nordstrom de Sitter black hole using the Poshl- Teller potential approximation. 
Nevertheless, in these papers the study of the Dirac QNMs were limited by using WKB or 
Poshl- Teller potential approximative methods. The reason to obstruct study of the Dirac 
QNMs is that for static black holes the standard wave-equations 

d 2 



possesses the special potentials 



h2 +^]Z ± = V ± Z ± II. I; 



v^x^±x±^ (1 . 2) 



which are the functions of a/A |2J| |2j| [2j| , where A is a function related to metric, say A = 
r 2 — 2Mr for the Schwarzschild black hole. So, it is very hard to calculate Dirac QNMs by using 
the known numerical methods (except WKB or Poshl- Teller approximative methods) since we 
have to expand the potentials as a series at boundary of the field, such as at the event horizon. 
However, Castello-Branco et al [2^ evaluated the Dirac QNMs of the Schwarzschild black hole 
by using the continued fraction approach in which they expanded the wave function as a series 
of y/A/r, and they found that the spacing for imaginary part is not —if AM, as it takes place 
for scalar, electrom agne tic, and gravitational perturbations, but —i/8M. 

Recently, we [3(3| |3jJ found that the wave functions and potentials of the Dirac field in the 
static spacetimes can be expressed as new forms, and the new wave functions are related to 
the standard wave functions Z± in Eq. (jl.lj) in a simple way j30|. Starting from the new wave 
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functions and potentials, wc [30] [31] showed that the Dirac QNMs of the Schwarzschild anti- 
de Sitter and Reissner-Norstrom anti-de Sitter black holes can easily be worked out by using 
Horowitz-Hubeny approach 0]. We have checked that the new potentials also present correct 
quasinormal frequencies for the scalar and gravitational perturbations in the Schwarzschild 
black hole. 

The main purpose of this paper is to study Dirac QNMs of the Schwarzschild black hole by 
using the new kinds of wave functions and potentials with the help of the continued fraction and 
Hill-determinant approaches. The results obtained here are different from Castello-Branco's 
one j29] (If taking M = 1, we can only obtain results n = 0, 2, ... in Ref. j^J. That is to say 
we do not find "specific modes".). We are sure that we do not miss any quasinormal frequency 
because the Hill-determinant approach gives all roots for n < 15. 

The organization of this paper is as follows. In Sec. 2 the decoupled Dirac equations and 
corresponding wave equations in the Schwarzschild spacetime are obtained by using Newman- 
Penrose formalism. In Sec. 3 the numerical approach to computing the Dirac QNMs is intro- 
duced. In Sec. 4 the numerical results for the Dirac QNMs in the Schwarzschild black hole are 
presented. The last section is devoted to a summary. 

II. DIRAC EQUATIONS IN THE SCHWARZSCHILD SPACETIME 

The Dirac equations [32j are 



where V bb' is covariant differentiation, P and Q are the two-component spinors representing 
the wave function, Pb> is the complex conjugate of Pb, and /x is the particle mass. In the 
Newman-Penrose formalism the equations become 



V2V B B'P B + i^Qs* = 0, 
V2V BB 'Q B + ¥Pb> = 0, 



(2.1) 



(D + e 
(A + /i 
(D + e 
(A + /Z 



p)P° + (5 + 7T- 
7 )P 1 + (5 + /? 
P)Q°' + (6 + 7T - 

7)Q 1 ' + (5 + ^ 



a)P l 
■t)P° 
a)Q v 

-f)Q l 



2- 1 /VQ 1 ', 

= -2- 1 / 2 z Ai P 1 , 
= 2- 1 /V^°- 



(2.2) 



For the Schwarzschild spacetime 




r 



(2.3) 
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where M represents the black hole mass, the null tetrad can be taken as 



P = (-, 1, 0, 0), 



with 



Then, if we set 



„ 1/ A 
n" = -(1, - -, 0, 0) 



"'" = 7Tr[°> °' ''sine 



A = r 2 - 2Mr. 



(2.4) 



(2.5) 



P 1 



-M_ 1/2 (r)5_ 1/2 (^)e- l ^-^ : 
M + i/ 2 (r)5 +1/2 (^) e - 1 ^-^, 
= M+i/2(r).S_ 1/2 (0)e- 1 ^-^, 



go = __R_ 1/2 ( r )5- +1/2 (0) e 



-i(u>t—mtp) 



where u; and m are the energy and angular momentum of the Dirac particles, after 
calculation Eq. ()2.2j) can be simplified as 



y/AVoR- 1/2 = (A + 2 / ur)v / AM + i/ 2 , 
APJ(VA]R +1/2 ) = (A - z/ir)R_ 1/2 , 



£l/2'S+l/2 



-xs. 



•1/2) 



AS 



+1/2- 



with 



9 zif n dA 
n <9r A A cfr ' 

= 1 1 

n <9r A A dr ' 

dm 
C ri = ™ + -t— 5 + n cot 0, 
ay sin 6* 



" n <9# sinfl 
K = r 2 u. 



+ n cot (9, 



We can eliminate S+1/2 (or S-x/2) from Eqs. (|2.9|) and ()2.10|) and obtain 



-1/2 

lw£d9 \ Sme d9 



-T,2 



m 



+ 2mscos8 + s 2 cos 2 9 



sin 2 9 



+ s + A s 



S„ = 0, 



(2.6) 
tedious 

(2.7) 
(2.8) 
(2.9) 
(2.10) 



(2.11) 



(2.12) 
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here and hereafter we take s = +1/2 for the case S. 



+1/2 



and s = —1/2 for £_ 



1/2 



>-l/2j> 



and v4+i/2 = A 2 — 2s and A-%/2 — A 2 . The angular equation ()2.12|) can be solved exactly and 
A s = (Z — s)(l + s + 1), where / is the quantum number characterizing the angular distribution. 
So, for both cases s = +1/2 and s = — 1/2 we have 

A 2 =(/ + £f. (2.13) 

We will focus our attention on the massless case in this paper. Therefore, we can eliminate 
M_i/2 (or \/AR +1 / 2 ) from Eqs. ()2.7|) and ()2.8|) to obtain a radial decoupled Dirac equation for 
AR + i/2 (or M_i/ 2 ). Then, Introducing an usual tortoise coordinate 

2 



r 



-dr 



(2.14) 



and resolving the equation in the form 



A -s/2 



(2.15) 



we obtain the wave equation 



dri 



(u 2 - V,)* t 



(2.16) 



where 



Vs 



A d 

4r 2 dr 



dr 



+ 



s 2 r 4 



d 
dr 



+ isur A 



dr 



+ 



A 2 A 



(2.17) 



Although the form of the potential is more complicated than standard wave-equations 



d 2 

dri 



+ 0J' 



A 2 ^ ± A^^ |25j|2Jj|28|, we will see that 



Z± = V±Z± with the potentials V± 
we can easily work out the Dirac quasinormal frequencies of the Schwarzschild black hole from 
Eqs. (J2.16j) and (J2.17j) by using continued fraction and Hill-determinant approaches. 



III. NUMERICAL APPROACHES TO COMPUTING DIRAC QUASINORMAL 
MODES 

For the Schwarzschild black hole, the QNMs are defined to be modes with purely ingoing 
waves at the event horizon and purely outgoing waves at infinity [7|. Then, the boundary 
conditions on wave function at the horizon (r = r+) and infinity (r — > +oo) can be 
expressed as 

{ y +) (3.1) 
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where k is the surface gravity of the black hole. 

A solution to Eq. (j2.16|) that has the desired behaviour at the boundary can be written in 
the form 



m=0 



r — Ta 



If we take r + = 1, the sequence of the expansion coefficient {a m : m 
by a three-term recurrence relation staring with clq = 1: 

a ai + p a = 0, 

a m a m+ i + (3 m a m + lma m -i = 0, m = 1,2, 



(3.2) 

1,2,....} is determined 



(3.3) 



The recurrence coefficient a m , j3 m , and 7 m are given in terms of m and the black hole parameters 
by 

= m 2 + (C + l)m + C , 



/3m 



-2m 2 + (Ci + 2)m + C 3 , 

m 2 + (C 2 - 3)m + C 4 - C 2 + 2, 



(3.4) 



with 



1 



Co 

Ci = -4 + 8zcj, 
C 2 = s + 3 — 4icj, 
C 3 = 8c; 2 - s - 1 + Uu - A 2 , 
4c; 2 



(3.5) 



C 4 = s + 1 - Auj' z - (2s + 4)iu. 

The three terms recurrence relation ()3.3|) can be used to find quasinormal frequencies by 
using both the continued fraction and Hill-determinant methods which are introduced in the 
following. 

[The continued fraction approach]: The ratio of the series coefficients {a m : m = 
1, 2, ...} is finite and can be determined in two ways: 

a m +i 1m a-m-i a m -27m-l a m -37m-2 



(X-m /3m— 1 /3m— 2 /3m—3 
—lm+1 Oim+llm+2 Ol m +2lm+3 



/3m 
O^m 



(3.6) 



/3 m +i~ /3 m +2~ (3 m +3~ 
Using Eqs. (|3.3|) and (|3.6|) we can obtain an equality between two continued fractions, one of 
infinite length, and the other finite: 



f3„ 



Olm-llm a m -2lm-l «07l 



m— 1" 



0; 



m-2" 



00 



Oimlm+1 Oim+llm+2 Ctm+2lm+3 



0m+l~ 0m+2~ 0; 



m+3~ 



(m = 1,2...). (3. 7) 
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This leads to a simple method to find quasinormal frequencies-define a function which returns 
the value of the continued fraction for an initial guess at the frequency. Then, use a root finding 
routine to find the zeros of this function in the complex plane. The frequency for which this 
happens is a quasinormal frequency. The nth quasinormal frequency is usually found to be 
numerically the most stable root of the nth inversion 34]. 

The success of the root-finding technique depends greatly on the initial guess given for the 
root. However, which initial approximation we shall take is a difficult work, especially at the 
beginning. Fortunately, we can easily find the low overtone quasinormal frequencies by means 
of the Hill-determinant method. 



[The Hill-determinant methodl: The condition that nontrivial solutions of Eq. (13.31) 

nn 

exist is given by the requirement that the Hill-determinant vanish 35][36J. The determinant 
can be easily seen in this case to be of the form 

A) "o 
7i fa ai 

72 A? "2 



D 



If D m stand for the (m + 1) x (m 
in terms of a difference relation 



7m- 1 /3m— 1 «m-l 
Tin Pm 



(3.8) 



1) determinant of D, then the D m can be concisely expressed 



(3.9) 



These equations are polynomials in u and can be solved at each order. However, the further 
requirement that the waves be purely outgoing at infinity and obey the stability conditions 
restrict the frequencies to complex plane with negative imaginary parts which are just the 
quasinormal frequencies we look for. 

The method is conceptually simpler than the continued fraction method of Leaver. However, 
the method looses accuracy for the higher overtones since one must evaluate the determinant 
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at values of order N (for giving N, one require solving polynomials of the order of 2N + 2 in 
uj) very much higher than the overtone number n of the quasinormal modes. Also, the method 
does not have any control over the stability of the root search algorithm. 

We have checked that both the methods present the same results for A = 5 to A = 10 and 
n = to n = 15. 



IV. DIRAC QUASINORMAL MODES OF THE SCHWARZSCHILD BLACK HOLE 

In this section we represent the numerical results of the Dirac quasinormal frequencies 
of the Schwarzschild black hole obtained by using the numerical approaches just outlined in 
the previous section. The results will be organized into two subsections: the fundamental 
quasinormal modes and high overtones. 

A. Fundamental quasinormal modes 

The fundamental quasinormal frequencies (n=0) for A = 1 to A = 40 are listed in the table 
(JTJ). Fig. (JTJ) shows that Au = — u\ as a function of A. From the table and figure we 
know that the fundamental quasinormal frequencies become evenly spaced for large A and the 
spacing is given by 

Au = 0.38490 - O.OOOOOi (4.1) 



B. High overtones 

The Dirac quasinormal frequencies of the Schwarzschild black hole for A = 1 to A = 5 and 
n = 1 to n = 20 are given by table (JT1J) and the first 21 Dirac quasi- normal frequencies for 
A = 1 to A = 10 are described by Fig. (J2J). 

The dependence of the quasinormal frequencies on A for the Schwarzschild black hole is 
described by Figs. © and (jlj). We learn from Figs. Q and (jlj) that A has the surprising effect 
of increasing real part of the quasinormal frequencies, but it almost does not affect imaginary 
part, especially for the lower overtones. 

From table (JTI]) and Fig. (JHJ) we know that the Dirac QNMs for A = 1, 2 demonstrate the 
following asymptotic behavior 

Re(u n+ i) — Re{oj n ) ^ as n — > oo, (4.2) 

i 

Im(u n+1 ) - Im(u n ) « -7T7 as n -> oo. (4.3) 



TABLE I: The fundamental Dirac quasinormal frequencies of the Schwarzschild black hole for A = 1 
to A = 40. 



A 


UJ 


A 


UJ 


A 


UJ 


A 


UJ 


1 


0.365926-0. 193965i 


11 


4.23212-0 


192462i 


21 


8.08197-0. 192453i 


31 


11.93128- 





192452i 


2 


0.760074-0. 19281i 


12 


4.61717-0 


192462i 


22 


8.46691-0. 192453i 


32 


12.31620- 





19245H 


3 


1.14819-0.1926H 


13 


5 00219-0 


192459i 


23 


8.85185-0. 192453i 


33 


12.70112- 


o 


19245H 


4 


1.53471-0.19254i 


14 


5 38720-0 


192457i 


24 


9.23679-0. 192453i 


34 


13 08604- 


o 


192451i 


5 


1.92059-0. 192507i 


15 


5.77220-0 


192456i 


25 


9.62172-0. 192452i 


35 


13.47096- 





192451i 


6 


2.30614-0.192490i 


16 


6.15718-0 


192456i 


26 


10.00667-0. 192452i 


36 


13.85588- 





192451i 


7 


2.69150-0.192479i 


17 


6.54215-0 


192455i 


27 


10.39160-0. 192452i 


37 


14.24080- 





192451i 


8 


3.07675-0. 192473i 


18 


6.92711-0 


192455i 


28 


10.77653-0. 192452i 


38 


14.62572- 





192451i 


9 


3.46192-0. 192468i 


19 


7.31207-0 


192454i 


29 


11. 16144-0. 192452i 


39 


15.01064- 





192451i 


10 


3.84704-0. 192464i 


20 


7.69702-0 


192454i 


30 


11.54636-0. 192452i 


40 


15.39554- 





192451i 



10 20 ^30 40 



10 20 30 40 

FIG. 1: The spacing Auj = uj\+i — UJ\ as the functions of A for the fundamental quasinormal 
frequencies. The top figure is drawn for AujR e = Re(uj\+i) — Re{uj\) which shows that the spacing of 
the real part is 0.38490 for large A. The bottom one is for Aujj m = Im(uo\ + \) — Im(uj\) which shows 
that the spacing of the imaginary part becomes zero for large A. 
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TABLE II: Dirac quasinormal frequencies of the Schwarzschild black hole for A = 1 to A = 5 and 
n = 1 to n = 20. 



n 


uj (A = 1) 


uj (A = 2) 


uj (A = 3) 


uj (A = 4) 




uj (A 


= 5) 


1 


0.295644-0.633857i 


0.711666-0.594995i 


1.11403-0.58543i 


1.50860-0.581936i 


1.89952-0.5802961 


2 


0.240740-1. 128450i 


0.638523-1.0369H 


1.05321-0.999391 


1.45954-0.98382i 


1.85899-0.9762311 


3 


0.208512-1.63397i 


0.569867-1.51491 


0.97949-1.44176i 


1.39383-1.40459i 


1 


80226-1 


38504i 


4 


0.187528-2. 14013i 


0. 516459-2. 01144i 


0.906879-1.91029i 


1. 3199-1. 8471H 


1 


7S409-1 


81 009i 


5 


0. 17248-2. 64556i 


0.476221-2. 51504i 


0.843213-2. 3969i 


1. 24587-2. 3100i 


1 


66006-2 


25269i 


6 


0.160962-3.15021i 


0.445224-3.0209i 


0.790087-2.89386i 


1.17721-2.78892i 


i 

± 




71 1 90i 


7 


0.151745-3.654191 


0.420575-3.527251 


0.746214-3.39616i 


1.11627-3.27894i 


I 


51418-3 


1851 3i 


8 


0. 144128-4. 15764i 


0.400401-4.03343i 


0.70972-3.90102i 


1.06324-3.77605i 


I 


AARAQ-1 


6691 6i 


9 


0.13768-4.66067i 


0.383487-4.53924i 


0.678957-4.4070H 


1.01733-4.27745i 


I 


38920-4 


1 6097i 


10 


0.132118-5. 16336i 


0.369026-5.0446i 


0.65265-4.9134H 


0.97747-4.78135i 


1 


33621-4 


658151 


11 


0.12725-5.66576i 


0.356461-5.54954i 


0.629847-5.41982i 


0.942637-5.286661 


1 


28901-5 


1 5892i 


12 


0.122937-6.16792i 


0.345399-6.05409i 


0.60984-5.92607i 


0.911955-5.79270i 


1 


24690-5 


66202i 


13 


0.119077-6.66989i 


0.335551-6.55828i 


0.59210-6.43207i 


0.884707-6.29908i 


1 


20921-6 


166621 


14 


0.115592-7.17170i 


0.326701-7.06215i 


0.57622-6.93778i 


0.860321-6.80556i 


1 


17532-6 


672151 


15 


0.112424-7.67335i 


0.318685-7.565741 


0.56189-7.44319i 


0.838336-7.311971 


1 


14467-7 


17822i 


16 


0.109525-8.174891 


0.311372-8.06907i 


0.548867-7.9483H 


0.818386-7.818251 


1 


11682-7 


684571 


17 


0.106858-8.676311 


0.304662-8.572181 


0.536959-8.453161 


0.800173-8.324351 


1 


09138-8 


191051 


18 


0. 104393-9. 17764i 


0.298471-9.07509i 


0.526008-8.95774i 


0.783457-8.83023i 


1 


06803-8 


69753i 


19 


0.102104-9.67888i 


0.292733-9.57782i 


0.515889-9.46207i 


0.768038-9.33589i 


1 


04651-9 


20394i 


20 


0.099979-10.1800i 


0.287392-10.0804i 


0.506497-9.966181 


0.753754-9.841321 


1 


02661-9 


710271 



It is shown that the spacing for imaginary part is — i/AM which is same as that of the scalar, 
electromagnetic, and gravitational perturbations. 



V. SUMMARY 

The wave equations for the Dirac fields in the Schwarzschild black hole spacetime are ob- 
tained by means of the Newman- Penrose formulism. Then, the quasinormal frequencies cor- 
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FIG. 2: First 21 Dirac quasi-normal frequencies of the Schwarzschild black hole for A = 1 to A = 10. 
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FIG. 3: Dependence of the real parts of the quasinormal frequencies on A for the Schwarzschild black 
hole, which shows that A has the surprising effect of increasing Re(u). 

responding to the Dirac field perturbation in the black hole spacetime are evaluated by using 
continued fraction and Hill-determinant approaches. Three interesting results obtained in this 
paper are listed in the following: (i) The fundamental quasinormal frequencies become evenly 
spaced for large A and the spacing is given by uj x +i — oo\ = 0.38490 — O.OOOOOi. (ii) A has the 
surprising effect of increasing real part of the quasinormal frequencies, but it almost does not 
affect imaginary part, especially for the low overtones, (iii) The quasinormal frequencies also 
become evenly spaced for large n and the spacing for imaginary part is —i/AM, which is same 
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FIG. 4: Dependence of the imaginary parts of the quasinormal frequencies on A for the Schwarzschild 
black hole, which shows that A almost does not affect imaginary part, Im(uj), especially for the lower 
overtones. 
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FIG. 5: Graphs of Aw = uj n +i — w n versus overtone number n for A = 1, 2. The top figure is drawn 
for AujRe = Re(uj n+ i) — Re(uj n ) which shows that Au)R e tends to zero for large n. The bottom one is 
for Au!j m = Im(uj n+ i) — Im(uj n ) which shows that Aw; m « — 437 f° r large n. 
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as that of the scalar, electromagnetic, and gravitational perturbations. 

In future papers in this series the new wave functions and potentials will be applied to study 
of Dirac QNMs in other black hole spacetimes, such as the charged Reissner-Norstrom black 
hole and the rotating Kerr black hole. 
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